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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2828 

EFFECT OF A FINITE TRAILING-EDGE THICKNESS ON THE DRAG 
OF RECTANGULAR AND DEITA WINGS AT SUPERSONIC SPEEDS 
By E. B. Klunker and Conrad Rennemann, Jr. 

SUMMARY 


The effect of a finite trailing-edge thickness on the pressure 
drag of rectangular and delta wings with truncated diamond-shaped air- 
foil sections with a given thickness ratio is studied for supersonic 
Mach numbers, linearized theory being used to evaluate the surface 
pressures. In order to facilitate comparison with wings having sharp 
trailing edges, the position of maximum thickness and base height are 
determined for least pressure drag as functions of a base-pressure 
parameter. Comparison is then made between the drag of these wings and 
similar wings with a sharp trailing edge for various aspect ratios and 
thickness ratios as a function of stream Mach number. The calculations 
of the drag characteristics for these wings show that significant drag 
reductions are possible under some conditions at high supersonic speeds. 
These drag reductions are relatively independent of aspect ratio for 
the rectangular wings but depend considerably on aspect ratio for the 
delta wings; the smaller aspect ratios show the larger drag reductions. 
Calculations of the spanwise distribution of drag are included to com- 
pare further the effect of a base on the drag for different aspect 
ratios . 


INTRODUCTION 


A great deal of attention has been focused on the problem of deter- 
mining supersonic profiles of min im um pressure drag. In reference 1, 
for example, profiles of minimum drag are determined by linearized theory 
for several auxiliary conditions, and in reference 2 an analysis is made 
for two structural conditions by using a nonlinear pressure relation. 

An interesting feature of these and similar investigations is that for 
certain conditions the profiles of minimum pressure drag have blunt 
trailing edges. The use of a finite trailing-edge thickness results in 
a reduction in the pressure drag of the forward surfaces of the airfoil 
(that is, the surface excluding the base) since the average absolute 
value of the slope of the airfoil surface is diminished. At high super- 
sonic speeds, the reduction in pressure drag of the airfoil surface may 
exceed the additional base drag incurred for some conditions. 
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Consequently, an airfoil with a finite trailing-edge thickness may have 
less drag than a wing with a sharp trailing edge which satisfies the 
same structural conditions. 

References 1 and 2 are concerned only with the effects of a finite 
trailing-edge thickness on two-dimensional airfoils. Of more practical 
interest is the effect of a "base on wings of finite aspect ratio, par- 
ticularly low- aspect -ratio wings, since it appears that these wings will 
he used at high supersonic flight speeds. An investigation of the drag 
of delta and rectangular wings is made in the present paper by using 
linearized theory to evaluate the surface pressures in order to deter- 
mine the effect of a finite trailing-edge thickness on the drag of 
supersonic wings. In order to simplify the analysis, a truncated 
diamond- shaped airfoil - the shape for minimum drag for a given thick- 
ness ratio for two-dimensional airfoils - is considered. The location 
of maximum thickness and the trailing-edge thickness are determined for 
least drag in order to reduce the number of parameters involved and to 
indicate the drag reductions possible. These optimum wings, however, 
are not wings of min imum drag since the wing plan form and the form of 
the profile are specified and only its parameters are determined to give 
minimum drag. The drag coefficients for these wings are compared with 
those for similar wings with a sharp trailing edge for several aspect 
ratios and thickness ratios for the Mach number range from 2 to 8. 


SYMBOLS 


c 



1 - r 
A 


chord 

root chord 
thickness ratio 
wing span 

ratio of base height to chord or ratio of base area to 
wing area 

position of maximum thickness, measured from leading edge, 
fraction of chord 

sweep angle of lea din g edge 

wing slope in stream direction 
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Xp slope of front surface in plane y = Constant 

/ 

Xq slope of rear surface in plane y = Constant 

M stream Mach number 

p Mach angle, sin -1 ^ 

P =\ju 2 - 1 

tan A _ _4_ 

P " PA 


U stream velocity 

w vertical disturbance velocity component 

0 disturbance velocity potential 


% 


x,y,z 


S 

Sb 

A 

P 


Cartesian coordinates 

0 

wing area 
base area 
aspect ratio 
pressure coefficient 

base pressure coefficient at any spanwise station 



average base pressure coefficient 
vacuum pressure coefficient 
base-pressure parameter 
wing pressure-drag coefficient 

pressure-drag coefficient of a wing with a sharp trailing 
edge 

section pressure-drag coefficient 
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c d section pressure-drag coefficient of a wing with a sharp 

trailing edge 


Subscripts : 

l wing of least drag 

0.5 wing with maximum thickness location at midchord 


AMLYSIS 

I 


The pressure drag of a wing is the drag due to the normal forces 
acting over its surfaces. The pressure drag of a wing with a finite 
trailing-edge thickness consists of a base drag and a pressure drag of 
the wing surface excluding the base. For supersonic flow, these two 
contributions to the pressure drag may be considered independently since 
disturbances from the trailing edge are not propagated upstream except 
to a small extent in the boundary layer. The normal pressures over the 
wing surface (excluding the base or trailing-edge surface) are deter- 
mined mainly by the flow field outside the boundary layer where the 
effects of viscosity are negligible, whereas the base pressures are 
related directly to the viscous wake. In general, the base pressures 
are determined from experiment, whereas the surface pressures may be 
determined from inviscid-flow theory. 

The surface pressures due to a thickness distribution may be calcu- 
lated by linearized theory by using the concept of source distributions. 
Puckett (ref. 3) has shown that for a thin symmetric wing in a super- 
sonic flow the required source strength at any point on the wing is pro- 
portional to the vertical component of the disturbance velocity ( slope 
of the wing in the stream direction) at that point. Thus the determi- 
nation of the disturbance velocity components, and consequently, the 
pressure, requires only an integration of the sources over the forward 
Mach cone. 

# 

The drag of a delta wing with a double-wedge airfoil section was 
calculated in reference 3 . The drag of a delta wing with a truncated 
diamond-shaped airfoil section may be calculated by slightly modifying 
the equations of reference 3- Since the calculation of the drag of a 
rectangular wing is the same in principle as that of the delta wing, 
many of the details of the analysis are omitted. These may be found in 
reference 3 and elsewhere. 

Rectangular wing.- For supersonic flow, the base drag may be con- 
sidered separately from the remainder of the pressure drag. Hence, the 
pressure -drag coefficient Gp of a thin symmetric wing whose airfoil 
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sections are all similar can be written as 



1 1 p »' as 


where P is the surface pressure coefficient, P^,' is the base pressure 
coefficient at any spanwise station, X is the wing slope in the stream 
direction, S is the wing area, is the base area, and ds is an 

element of area. If the average spanwise value of the base pressure 
coefficient is denoted by F^, the wing pressure drag may be written as 



( 1 ) 


and the quantity h is the ratio of the base height to the chord, or, 
since all the airfoil sections for the wing are similar, it is the ratio 
of base area to wing area. The pressure drag of: the wing surface 
excluding the base is considered in several parts for convenience. 

The supersonic flow field -past a rectangular wing is two-dimensional 
except in the region bounded by the wing tip and the tip Mach cone: 



■>- x, i 
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The surface pressure coefficient for the two-dimensional part of the 
flow field (region I, composed of the subregions 1^ and Ig of the 
sketch) is 

P = 2 £ 

0 


where 0 = yM^ - 1 and M is the stream Mach number. The drag coef- 
ficient of this two-dimensional-flow region then is 

^ ■ w f Sz x2i ° 


and for a double-wedge airfoil section with a finite trailing-edge 
thickness 




(1 - r) 



where 1 - r is the location of maximum thickness measured from the 
leading edge and A is the aspect ratio. 

The wing plan forms and airfoil sections for both the rectangular 
and delta wings are shown in figure 1. From this figure the slope b of 
the forward and rear surfaces are found as 


h. ~ 


t 

2(1 - r) 




where t is the thickness ratio and B = The drag coefficient for 
the two-dimensional region then becomes 


CD! = 


iiJ 

1 

1 - B(1 - r)(2 - B) 1 

, . (1 - B) 2 (2 - r) ] 

[3 

r(l - r) 0A 

t 

L 1 + - Jj 


(2) 


for BA > 2. The condition BA > 2 requires that the tip Mach cones do 
not intersect on the wing surface. 
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The pressure coefficient at a point (x,y) in the tip region (com- 
posed of the subregions Hi; Il2.> H 3 of the sketch ) } which is 

bounded by the wing tip and the tip Mach cone, is found from the poten- 
tial of the source distribution in the forward Mach cone. In refer- 
ence 3 the required source strength at any point is shown to equal the 
vertical component of velocity w at that point. The velocity poten- 
tial 0 at (x;y,0) in the tip Mach cone for a constant source strength 
is 




ny px-p(y-rj) 

0 (x,y,O) = -7 / drj / - 7= 

0 J ° \/(* - I ) 2 - P 2 (y - i )) 2 


x+Py 


JJT- f 


+P(y-n) 


dl 


\/(x - S) 2 - p 2 (y - tj ) 2 


and integrating with respect to S gives 

'X+ffy 


0(x,y,O) =-2 f P cosh" 1 - — r dq 

« J 0 P|y - n) 


The component of the disturbance velocity in the x-direction is found 
after some reduction as 


0 x (x,y,O) 




and from the linearized form of the Bernoulli equation the pressure 
coefficient is 


P 




(3) 


where U is the stream velocity and the airfoil slope is 

The pressure coefficient for two-dimensional flow and the pressure 
coefficient given by equation (3) can be used to calculate the drag of 
the tip regions. The tip drag is determined by considering sources of 
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strength \]_U distributed over the whole tip area and sources of 
strength ( - Xi)U distributed over the regions II 2 and H 3 . For 
the sources of strength Xq_U the pressure coefficient is given by- 
equation ( 3 ) and the contribution to the drag for both wing tips is 



For sources of strength (Xg - Xi)U, the pressure coefficient is the 
two-dimensional value in the region II 2 (since this region is outside 
the Mach cone emanating from the maximum-thickness line) and the pressure 
coefficient is given by equation ( 3 ), with x replaced by x - c(l - r), 
for the region II 3 . This contribution to the drag of both tips is 



2t 2 (l - B) [l - B(1 - r)] t 2 (l - B)[l - B(1 - r)] it - 1 


P 2 Ar 


P 2 A(1 - r) 


it 


(5) 


The total pressure drag of the rectangular wing then is 


C D “ + C D X + C D 2 “ P b h 


or, with h = Bt, 


C D ” C D + V 


B(B - 2) 


-2i- 

«PA t 


(PA > 2 ) 


( 6 ) 
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where 


C 


D 


1 

P r(l - r) 


is the drag of a rectangular wing with a sharp trailing edge. It is of 
interest to note that the drag of a rectangular wing with a finite 
trailing-edge thickness depends upon the aspect ratio; whereas the drag 
of a wing with a sharp trailing edge is independent of the aspect ratio. 

Delta wing .- The pressure drag of a delta wing with a double-wedge 
airfoil section has been calculated in reference 3 * The drag of a sim- 
ilar wing with a finite trailing-edge thickness can be found from the 
same relations provided the expression for the airfoil slope on the rear 
surface is suitably modified and the base drag is added. Since the 
effect of a finite trailing-edge thickness on the drag reduction is 
significant only at relatively high supersonic Mach numbers, only the 
condition of supersonic leading edges is presented. This corresponds 
to case 1 of reference 3 . 


The pressure drag of the delta wing C-^, excluding the base drag, 
is, from equations (4l), (42), and (43) of reference 3, 

C Di . G l( n,r) ♦ ?^£[f - <*(»,*>] + ^ r § ( 7 ) 


where 


G x (n,r) = 


1 - r 


1 + r 


\Jl - n c 


cos _1 n + 


yi - n^r^ ' 


lie . -1 
— + sin nr 


and n is the ratio of the tangent of the leading edge sweep angle A 

to the cotangent of the Mach angle, that is, n = - a fi — A . = it. The pres- 

P PA 

sure drag may be expressed as the sum of the drag contribution from 
equation (7) and the base drag. Hence 


C D = C D 1 ~ Pbk 
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where the average spanvise value of the base pressure coefficient P-^ 
for the delta wing is 


Pb = 2/ 1 V #- 

J 0 c r V 3 / 2 / 


where c is the local chord and c r is the root chord. By substi- 
tuting from equation (J) the drag coefficient becomes 


C D = C D [l - B(1 - r)] - ^ 


t 2 

B(1 - B) 

p 

r 


+ B 




where 


C-_ t2 
D “ R 


2 

cos“^n 

ir(l - r 2 ) 

\jl - n 2 


r(~ + 


sin 


•^r) 


( 8 ) 


is the drag of a delta wing with a sharp trailing edge and n < 1. 

Lea st -drag wings .- The equations for the drag of both the reotan- 

PPb 

gular and delta wings each contain four parameters: r, B, ", 

. t 

and PA (or n) . The effect of a finite t railing-edge thickness on the 

drag can be more clearly demonstrated by reducing the number of param- 
eters. If the position of maximum thickness and the base height are 
determined to give the least drag for a given thickness ratio, then the 

PP b 

quantities r and B may be expressed in terms of the parameters — 

and PA (or n) . The wings so determined, however, are not wings of 
minimum drag for a given thickness ratio in an exact sense since the 
wing plan form and the form of the profile are specified and only its 
parameters are determined to give minimum drag. The wings so determined 
are herein denoted as least -drag wings. 

The values of r and B required to minimize the drag are found 
from the simultaneous solution of the equations 


dr 


= . 0 


dC 


D 
dB 


= 0 


(9) 
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in terms of the quantities -- 


and pA (or h) . These values may 


then be substituted into equation (6) or (8) to yield the least drag 

6Ph 

coefficient for the particular wing- plan form as a function of — —a 

u 

and PA (or n) . 


DISCUSSION 


The wing profile geometry which minimizes the drag of a given wing 
plan form with plane surfaces is found as the simultaneous solution of 
equations ( 9 ) . It is significant that the least -drag conditions for the 
rectangular wing are similar to those for a two-dimensional airfoil. 

The two-dimensional mini mum -drag airfoil for a given thickness ratio, 
as determined by linearized theory, satisfies the .condition = -X^y 
further, the position of maximum thickness and the trailing-edge thick - 


pPb , 

ness depend only on the base-pressure parameter — — (see ref. 2, for 

X 

example) . For the rectangular wing the condition = ~^2 holds 

for the least -drag condition (see appendix) however, the position of 
maximum thickness and the trailing-edge thickness depend both upon the 
base -pressure parameter and the aspect -ratio parameter PA. For the 
delta wing the two slopes are not equal and the geometric parameters 

0 p b 


for least drag are functions of — 


and 


The results of the solution of equations (9) are presented in 
figures 2(a) and 2(b) for the rectangular wing and the delta wing with 
supersonic leading edges. In these figures the position of maximum 
thickness 1 - r and the ratio of trailing-edge thickness to wing 

PPk 

thickness B are given as functions of the parameters — ~ and PA 

(or n) . These values may be substituted into equations (6) and (8) to 
determine the least drag coefficients. In figures 3(a) and 3(b) the__ 
ratio of the least drag coefficients CL to the drag coefficient CL. 

is the drag coefficient of a similar wing with a double- 


( where C 


d 0.5 


wedge airfoil section and maximum-thickness location at midchord) is pre- 
sented as a function of for the rectangular and delta wings, 

U 


respectively. Figures 2(a) and 3(a) show that for the rectangular wing 
both the geometric parameters and the drag coefficients differ but little > 
(for moderate values of PA) from the values for the two-dimensional 
flow. Further, large pressure -drag reductions appear possible with the 
use of a finite trailing-edge thickness. Whether these large drag reduc- 
tions can actually be realized, however, depends upon the base pressures 
experienced . 
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In order to assess the value of wings with a finite trailing-edge 
thickness, the elimination of the "base pressure as a parameter is neces- 
sary. For this purpose, the base pressure coefficients presented 'in 
reference 2 have been utilized and are presented in figure 4 as a func- 
tion of Mach number. Vacuum pressure coefficients are also presented 
for comparison. These values of the base pressure are representative 
for turbulent boundary layers. 


CD7 


The ratio of the drag coefficients — 


for the base pressures 


of figure 4 are presented in figures 5(a) and 5(^) for the rectangular 
and delta wings, respectively, as a function of the stream Mach number. 

As in the case of the two-dimensional airfoil, the percentage drag reduc- 
tion possible with the use of a finite trailing-edge thickness is depend- 
ent to a large extent on the wing thickness ratio; the thicker wings show 
a substantially larger drag reduction at a given stream Mach number.-^* 

For example, the drag ratios for the rectangular wing at a Mach number 
of 6 are approximately 0.76 and O.98 for thickness ratios of 0.10 
and 0.06, respectively. For the rectangular wing, the effect of aspect 
ratio on the drag reduction is small compared with the effect of thick- 
ness ratio; for a given Mach number and thickness ratio the difference 

for aspect ratios of 1 and °° is generally less than 3 per- 


in 


- ^ 0.5 

cent. In contrast, the effect of changes in thickness ratio and aspect 
ratio are of the same order of magnitude for the delta wing and the 


difference between 


than 10 percent. 


^ 0-5 

The ratio 


ratio for both the rectangular and delta wings . 


for aspect -ratios of 1 and 4 may be more 

C D 1 

■ • ■ ■ decreases with decreasing aspect 

Cd o.5 


A clearer understanding of the effect of aspect ratio on the drag 
of the rectangular and delta wings is obtained by considering the span- 
wise distribution of section drag. The spanwise distribution of the 
quantity cc,j, where c is the section chord and c<j is the section 
drag coefficient, for these wings is presented in figure 6 for a stream 
Mach number of 5 and a thickness ratio of 0.10. For simplicity, the 
base pressure coefficient is taken as constant over the span. 


^The drag coefficient of the two-dimensional airfoil of minimum 
drag c dm1r| is a linear function of t, whereas the drag coefficient of 

the airfoil of fixed geometry cjj increases as t^ and the ratio • 
decreases with increasing values of the thickness ratio. C( ^ 
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The pressures "both inside and outside the Mach cone are unaffected 
by a change in aspect ratio for the rectangular wing with a specified 
airfoil section. Thus, changing the aspect ratio only changes the pro- 
portion of the two-dimensional region to the tip region and the spanwise 
distribution of drag changes similarly. Sihce the airfoil geometry for 
least drag depends slightly upon the aspedt ratio, the pressures, and, 
consequently, the distribution of drag for the wing of least drag vary 
somewhat with aspect ratio. Figure 6 shows that the drag distributions 
in the tip regions are similar for different aspect ratios. That is, 
the basic distribution of drag within the tip does not change as the 
aspect ratio is changed although the proportion of tip area to wing area 
changes with, a change in aspect ratio. Moreover, in a large part of the 
tip region, the difference between cl and c= is greater than in 

0.5 Cj) 7 

the two-dimensional region. Consequently, the ratio of = — — decreases 

Cd 0.5 

with decreasing aspect ratio since the tip region is a larger percentage 
of the total wing area for the lower aspect ratios. 


A change in aspect ratio for the delta wing necessarily alters the 
leading-edge sweep angle, and the pressures both inside and outside the 
Mach cone from the apex are changed. As shown in figure 6, a change in 
aspect ratio also alters the form of the spanwise distribution of drag 
for the airfoil of fixed geometry. For the delta wing of least drag, 
however, the distributions for the different aspect ratios tend to be 
similar; the section drag approaches a linear variation across the span. 


CONCLUDING REMARKS 


The calculations of the drag characteristics for rectangular and 
delta wings for a given thickness ratio show that significant drag 
reductions are possible under certain conditions at high supersonic 
speeds with the use of a finite trailing-edge thickness. The Mach num- 
ber at which a blunt-trailing-edge wing has less drag than a similar 
wing with a sharp trailing edge is dependent upon the wing plan form, 
aspect ratio, and thickness ratio. For. a given Mach number and aspect 
ratio, the thicker wings show a substantially larger drag reduction 
through the use of a finite trailing-edge thickness. The effect of 
aspect ratio on the drag reduction for the rectangular wings is small, 
whereas it is of the same order of magnitude as the effect of thickness 
for the delta wing. The ratio of the least drag coefficient of a wing 
with a blunt base to the drag coefficient of a similar wing with a sharp 
trailing edge decreases with decreasing aspect ratio. For the rectan- 
gular wing the spanwise distributions of drag in' the tip region are 
similar for all aspect ratios, although the proportion of tip area to 
wing area changes with changing aspect ratio. For the delta wing with 
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a fixed airfoil section, the form of the spanwise drag distribution 
is altered as the aspect ratio is changed, whereas for the least-drag 
delta wings the distributions for different aspect ratios tend to be 
similar. Although the analysis of this paper is confined to rectangu- 
lar and delta wings with a given thickness ratio, the results suggest 
that significant drag reductions my be possible by the use of a finite 
trailing-edge thickness with other plan forms. It my be expected that, 
similar to the result for the two-dimensional airfoil (KACA TN 2264), 
the relative drag reduction possible by the use of a finite trailing- 
edge thickness with other structural conditions (given area, torsional 
stiffness, and so forth) is greater than for a given thickness ratio. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 15, 1952. 
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APPENDIX 

LEAST-DRAG CONDITIONS FOR RECTANGULAR WING 


The airfoil of minimum drag for a given thickness ratio in a two- 
dimensional supersonic flow is a truncated double wedge . Within the 
scope of linearized theory, is then equal to that is, the 

slopes over the front and rear surfaces are equal (see ref. 2, for 
example) . It is interesting to note that this same condition also 
applies to the rectangular wing of least drag with truncated diamond- 
shaped airfoil sections. 

The drag of a rectangular wing with a truncated diamond -shaped 
airfoil section is (eq. (6)) 



t2 

B(B - 2) 

B 2 

e p b 

0 

, r 

JT0A 

t 


B 


where 



£_L_ 

0 r(l - r) 


The values of . r and B required to give least drag are found 
from the simultaneous solution of the equations 


d c D 
8 r 


= 0 



(Al) 


From the first of equations (Al) 



2 2 


(A2) 


r jt£A 
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and from the second of equations (Al) 


B = 


1 - 2r 
1 - r 


(A3) 


After eliminating B from equations (A2) and (A3), r is found as 


2 + 


r = 




-(-$) 


k 

JtpA 




m 


The slopes of the front and rear surfaces are Xq. = 


2(1 - r) 


and 


respectively. Substituting equation (A3) for B in 
the equation for X2 results in the relation 


.. t L 1 - 2r\ t , 

^ « "2rf" 1 - r j “~2(1 - r) 1 


Consequently, the condition for least drag of the rectangular wing with 
a truncated diamond-shaped airfoil Beet ion is Xq_ = -Xg. 

From equations (A2) and (A4) it is seen that r and B depend on 

f^b 

both the base pressure parameter and the aBpect ratio pa ram- 

"fc 

eter PA* Tills is in contrast with the two-dimensional case where r 

pPh 

and B depend only on - 
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Figure l.~ Wing geometry for rectangular and delta wings 
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(a) Rectangular wing. 

Figure 2.- Position of maximum thickness and trailing-edge thickness for 
minimum drag as a function of the base-pressuri pSameter 
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(b) Delta wing. 
Figure 2.- Concluded 


Baee-preaaure parameter, H ■ - 
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(a) Rectangular wing. 

Figure 3.- Ratio of mi nira i mi drag coefficient of a vdng with a blunt 
trailing edge to drag coefficient of a wing with a double-wedge 

R#ar»4’.*i rin n o o A /■**-» r\-f* 4* Vi /-> v.vws i 
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Base-pressure parameter, H 


(b) Delta wing. 
Figure 3*- Concluded. 
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Figure U.- Vacuum pressure coefficient and base pressure coefficient as 

functions of stream Mach n umb er. 
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Drag -coefficient ratio 
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Stream Haoh mmber, M 


(a) Reciaugulaj.- wing. 


Figure Ratio of minimum drag coefficient of a wing with a blunt 
tr ailing edge to drag coefficient of a wing with a double-wedge 
section as a function of the stream Mach number for base pressures 
presented in figure 4- 
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.(b) Delta wing. 
Figure 5.~ Concluded. 
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Distanoe along Bpan, 


Distanoe along Bpan, 


(a) Rectangular wing. (b) Delta wing. 

Figure 6.- Comparison of epanwise distribution of drag for rectangular and 
delta wings of least drag with rectangular and delta wings having sharp 
trailing edge and position of maximum thickness at midchord. 11 = 
t “ 0.10. 
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